THE TEMPERATURE-MASS RELATION IN MAGNETIZED GALAXY CLUSTERS 



Klaus Dolag 1 , August Evrard 2 , and Matthias Bartelmann 1 
'Max-Planck-Institut fur Astrophysik, P.O. Box 1317, D-85741 Garching, Germany 
2 Dept. of Physics, University of Michigan, Ann Arbor, MI 48109-1 120, USA 



o 
o 

(N 

c 

a 



ABSTRACT 

We use cosmological, magneto-hydrodynamic simulations of galaxy clusters to quantify the dynamical impor- 
tance of magnetic fields in these clusters. The set-up of initial magnetic field strengths at high redshifts is chosen 
such that observed Faraday-rotation measurements in low-redshift clusters are well reproduced in the simulations. 
We compute the radial profiles of the intracluster gas temperature and of the thermal and magnetic pressure in a set 
of clusters simulated in the framework of an Einstein-de Sitter and a low-density, spatially-flat CDM cosmological 
model. We find that, for a realistic range of initial magnetic field strengths, the temperature of the intracluster gas 
changes by less than w5%. 



1. Introduction 



y—{ Magnetic fields are common in galaxy clusters. They are in- 
ferred from Faraday-rotation measurements (Vallee et al. 1986, 
£N| 1987) in polarised background radio sources, and from diffuse 
CD radio haloes around galaxy clusters (Kronberg 1994). Depend- 
ed ing on several assumptions, most notably on the characteristic 
y—i length scale of field coherence, magnetic field strengths inferred 
CD from observations in typical galaxy cluster centres range from 
*~ 1 s» 0.1/jG up to some tens of /jG. 

CD Provided the magnetic field is tangled on sufficiently small 
^ <— | scales for it to be considered randomised, it contributes a non- 
O fhermal pressure p%, which is related to the thermal pressure p± 
I through 
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Conversely, magnetic fields ordered on scales comparable to the 
cluster scale can add anisotropic non-thermal pressure. 
■ The structure of intracluster magnetic fields is largely un- 
clear at present, and consequently the magnetic field strengths 
in clusters are equally uncertain. We investigate in this paper the 
amount of non-thermal pressure support due to realistic mag- 
netic fields in simulated galaxy clusters. 

Our approach is as follows (see Dolag et al. 1999). We set up 
gas-dynamical cosmological cluster simulations at high redshift 
and impose magnetic fields of variable strength and structure 
as summarised in Table ^ below. Since the origin of magnetic 
fields on cluster scales is unknown, we use two extreme types of 
seed fields, namely either completely homogeneous or chaotic 
initial magnetic field structures. The evolution of the intracluster 
fields is then computed during cluster collapse. The initial mean 
magnetic field strength (or energy density) is a free parameter 
in our simulations. We choose its range such that the simulated 
clusters at low redshift well reproduce the statistics of observed 
Faraday-rotation measurements taken from Kim et al. (1991). It 
turned out in earlier work (Dolag et al. 1999; hereafter DBL99) 
that (i) the initial field structure is irrelevant because the final 
field structure is completely determined by the dynamics of the 
cluster collapse, and (ii) initial field strengths of order 10~ 9 G at 
the initial redshift of the simulations (see Table |l]) are adequate 
to reproduce the measured Faraday-rotation statistics. 



2. Models 

We use a set of 90 numerically simulated galaxy clusters for in- 
vestigating the effect of magnetic fields on the temperature and 
pressure stratification in these clusters. Here, we summarise the 
numerical technique only briefly because details on the meth- 
ods applied, their implementation, and their performance in test 
problems have been published elsewhere (DBL99). 



2.1. Numerical method: GrapeMSPH 

We use the cosmological MHD code described in DBL99 to sim- 
ulate the formation of magnetised galaxy clusters from an initial 
density perturbation field. The code combines the merely gravi- 
tational interaction of a dark-matter component with the hydro- 
dynamics of a gaseous component. The gravitational interac- 
tion of the particles is evaluated on GRAPE boards (cf. Sugi- 
moto et al. 1990), while the gas dynamics is computed in the 
SPH approximation (Lucy 1977; Monaghan 1992). The original 
"GrapeSPH" code (Steinmetz 1996) was extended by adding the 
magneto-hydrodynamic equations to trace the evolution of the 
magnetic fields. Because of the assumed ideal electric conduc- 
tivity, the fields are frozen into the gas flow. The back-reaction 
of the magnetic field on the gas via the Lorentz force is included. 
The numerical viscosity required to capture shocks in SPH is 
chosen such that angular-momentum transport in presence of 
shear flows is carefully controlled. The SPH kernel width is 
automatically adapted to the local number density of SPH parti- 
cles, which results in an adaptive spatial resolution of the code. 

Extensive tests of the code were performed and described in 
DBL99. The code succeeds in solving the co-planar MHD Rie- 
mann problem posed by Brio & Wu (1988). Although the simu- 
lated magnetic field is not strictly divergence-free, V B is always 
negligible compared to the magnetic field divided by a typical 
length scale of B (see DBL99 for details). The code also assumes 
the intracluster medium to be an ideal gas with an adiabatic in- 
dex of y = 5/3, and neglects cooling. The surroundings of the 
clusters are important because their tidal fields affect the overall 
cluster structure and the merger history of the clusters. There- 
fore, the cluster simulation volumes are surrounded by a layer 
of boundary particles whose purpose it is to accurately represent 
the tidal fields of the cluster neighbourhood. Here we specially 
focus on the dynamical influence of the magnetic field in galaxy 
clusters, e.g. on the pressure support and the temperature in the 
clusters. 



2.2. Initial conditions 

We set up cosmological initial conditions for two CDM- 
dominated universes (EdS and FlatLow) as specified in Table |l} 
For each cosmology, we calculate ten different realisations of 
the initial density-fluctuation field at redshift z\m, which result 
in clusters of different final masses and dynamical states at red- 
shift z — 0. We simulate each of these clusters with up to five 
different initial magnetic field configurations, listed in Table g, 
yielding a total of 90 cluster models. 



Table 1 . — Parameters of the EdS and the FlatLow models. 



model 


H() Q.Q A 0"8 /baiyon £ini 


EdS 

FlatLow 


0.5 1.0 0.0 1.2 5% 15 
0.7 0.3 0.7 1.05 10% 20 



In lack of any detailed knowledge on the origin of primordial 
seed fields, we explore two extreme cases of initial field con- 
figurations. In one case ("homogeneous"), we assume that the 
field is initially constant throughout the cluster volume. In the 
other case ("chaotic"), we let the initial field orientation vary 
randomly from place to place, subject only to the condition that 
V B = 0. The initial field strengths in both cases are determined 
by setting the mean field energy densities. Table || summarises 
the initial field set-ups and the mean field strengths in the cluster 
cores. The initial field strengths are of order 10~ 9 G, the final 
field strengths of order /jG. 



Table 2. — Initial magnetic fields (column 2) and the resulting 
final magnetic field strengths in the simulated clusters (column 3 
and 4 for EdS and FlatLow models, respectively). The final val- 
ues are an average of the magnetic field over the central region 
(within a radius of 350kpc) and over all ten clusters for each 
cosmology. 



initial field 


Bini 


(-Bfinal)core 


\HatLow 
\t> final /core 


no 


0.0G 






low (chaotic) 


0.2 x 10~ 9 G 


0.4 ^G 




low 


0.2 x 10~ 9 G 


0.4^ 


0.3 nG 


medium 


1.0 x 10~ 9 G 


1.1 fiG 


Q.%HG 


high 


5.0 x 10~ 9 G 


2.5 <uG 


2.0^G 



2.3. Evolution of intra-cluster fields 

At low redshifts, our simulated clusters reach temperatures of 
kT ~ lOkeV, in agreement with observed clusters of similar 
mass. The gas densities in their cores range within 10~ 4 — 
10~ 3 cm~ 3 , as inferred from X-ray luminosities and tempera- 
tures of observed clusters. 

In agreement with our earlier paper (DBL99), we find that 
the structure of the intra-cluster magnetic fields at low redshifts 
is entirely independent on the structure of the initial magnetic 
field. For homogeneous and chaotic initial field set-ups, the 
mean strength and the coherence length are the same. Initial 
fields of 10~ 9 G at Zini turn into final fields of /jG strength at red- 
shift zero, and the typical coherence length of the fields is of or- 
der 50kpc//i. In comparison to expectations based on magnetic 



flux conservation in spherically collapsing spheres, the amplifi- 
cation of the field is one order of magnitude higher because of 
shear flows in the intracluster plasma. 

During cluster collapse, the magnetic field is highly com- 
pressed, twisted, tangled and sheared. Any information on its 
initial structure is lost in that process. This is confirmed by 
Fig. [j], where it is shown that the magnetic pressure support does 
not depend on the initial field structure. We can therefore restrict 
our study to the homogeneous initial field set-up, which is nu- 
merically easier to handle. 

Due to the violent motion of the intracluster plasma during 
cluster collapse, the final magnetic field strength depends on the 
dynamics of the cluster formation. Starting with identical initial 
magnetic field set-ups, final field strengths differ by a factor of 
4 — 5, depending essentially on cluster mass. The mean mag- 
netic field strength drops by about two orders of magnitude be- 
tween the cluster centre and the virial radius. A detailed study of 
the evolution and the structure of intracluster magnetic fields in 
different cosmologies will be described in a forthcoming paper. 

A direct way to compare simulated magnetic fields in galaxy 
clusters with observations is provided by Faraday-rotation mea- 
surements. We created samples of synthetic rotation measures 
by integrating the Faraday rotation along 10 4 light rays ran- 
domly shot through each of a sub-set of our cluster simulations. 
As in DBL99, the cumulative rotation-measure distributions ob- 
tained in simulations with either homogeneous or chaotic initial 
field set-up were statistically identical, confirming once again 
that the structure of the initial field has no observable conse- 
quences for the structure of the final field processed by cluster 
collapse. 

We further performed two types of comparison between syn- 
thetic and observed Faraday-rotation data. In the first, we com- 
pared the radial distribution of rotation measures seen in the 
Coma cluster (Kim et al. 1990) with that obtained from one of 
our simulated clusters with mass similar to Coma. In the second, 
the sample of rotation measures compiled by Kim et al. (1991) 
from observations of many Abell clusters was compared to syn- 
thetic rotation measures obtained from our set of simulated clus- 
ters. In both cases, synthetic rotation-measure distributions were 
obtained along 10 4 light rays randomly shot through the simula- 
tions. 

Both comparisons showed that the medium initial mag- 
netic field strength could reproduce the observed statistics very 
well. In both cases, the simulated and the observed rotation- 
measure distributions could not significantly be distinguished in 
a Kolmogorov-Smirnov test. However, the scatter in both mea- 
surements and simulations is large. This is due to the relatively 
low number of simulated clusters per set, and to the selection 
criteria of galaxy clusters involved in the simulations and the 
observations. We therefore decided to allow for a range of ini- 
tial magnetic field strengths which results in a range of synthetic 
Faraday-rotation distributions which are still well in agreement 
with the observed data. Further detail can be found in DBL99 
and Dolag (2000). Comparing with more recent data by Clarke 
et al. (1999), we find identical or even improved agreement. A 
study of clusters in different cosmologies, using newer observa- 
tional data, is under way. 

Choosing a realistic model to distribute relativistic electrons 
in the simulated clusters, the simulations are also able to match 
the properties of observed radio-haloes very well, including the 
observed very steep relation between cluster temperature and ra- 
dio luminosity (see Dolag & Ensslin 2000). 




Fig. 1 . — Shown here is the ratio between the magnetic and the thermal pressure in shells around the cluster centre as a function of radius in units 
of the virial radius ^200 (heavy lines). Different line types distinguish the different initial magnetic field strengths as indicated in the plot. The 
two dotted lines represent the low field with homogeneous and with chaotic initial field configurations. Both lead to the same amount of pressure 
support in the final clusters. The curves show the average across the entire simulated cluster sample. The thin lines indicate the rms scatter within 
the sample. Near cluster centres, the ratio between magnetic and thermal pressure is consistent with the isotropic estimate of eq. (jlj). The left panel 
represents the realisations for the EdS cosmology, the right panel for the FlatLow cosmology. 



2.4. Quantities calculated 

The correlation length of final magnetic field configurations near 
cluster centres is of order 50kpc (see DBL99). Field orien- 
tations in regions large compared to the correlation length can 
therefore be considered as randomised. For instance, there is no 
preferred direction in the final magnetic pressure distribution. 
Thus, we can safely average thermodynamic quantities like the 
gas pressure and gas temperature in spherical shells, provided 
these shells are sufficiently thick compared to the field correla- 
tion length. 

The radial distances to the cluster centres are scaled by the 
virial radius r2oo- For simulated clusters, r2oo is defined as the 
radius enclosing a mean density p(r2oo) of 200 times the critical 
density p cr , 
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The temperature profiles give mass-weighted temperatures av- 
eraged in shells, divided by the mass-weighted temperatures in 
the shell at the radius rioo of the corresponding non-magnetised 
cluster. 



3. Results 

During cluster formation and evolution, the simulated intraclus- 
ter magnetic field strengths are amplified by compression, tan- 
gling and bending of field lines, and shear flows in the gas. The 
magnetic field structure is also affected by merging, accretion 
of sub-clumps, and consequent relaxation processes. In com- 
bination, these processes amplify the initial magnetic fields in 
the simulation by factors of order « 1000. Final fields reach 



strengths of a few /jG in and near cluster centres, in good agree- 
ment with field strengths inferred from observations. The grow- 
ing magnetic fields affect the balance between the gravitational 
force and the total (magnetic plus thermal) pressure during clus- 
ter formation. 

3.1. Radial structure 

On the whole, the fairly isotropic magnetic pressure has the ten- 
dency to counteract the gas collapse, and therefore acts in a sim- 
ilar manner as the thermal pressure. This implies that the gas 
stratification can achieve hydrostatic equilibrium at lower tem- 
perature, because lower thermal pressure is required in presence 
of the additional pressure support due to the magnetic field. Fig- 
ure [j] shows, averaged in shells, the ratio 
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between the magnetic pressure ps and the thermal pressure p± 
as a function of the radius R. Due to compression and turbu- 
lence in the gas flow, the magnetic field is highest in and near 
the centres of the simulated clusters, where the thermal pres- 
sure is also highest. As Fig. [l] shows, the relative contribution 
of magnetic fields to the pressure support appears to peak some- 
where between the cluster core radii and the virial radii r2oo- 
This feature, however, is not very pronounced, so that the ratio 
between magnetic and thermal pressure support is almost con- 
stant across a large range of cluster-centric radii. Due to the 
smaller final masses of the clusters in the FlatLow cosmology, 
the magnetic pressure support is also smaller in these clusters. 
Nevertheless, these less magnetized clusters are able to match 
the observed Faraday-rotation. This is because the value of Hq 




Fig. 2. — This figure shows the radial temperature profile, averaged across the simulated cluster sample. The temperature is scaled by the mass- 
weighted temperature of the non-magnetised cluster at a radius of r2oo, and the radius is scaled by ^200- Different line types distinguish different 
initial magnetic field strengths, as indicated in the plot. The temperature change with increasing magnetic field strength is clearly visible. The 
temperature decrease near cluster centres is compatible with the isotropic estimate of eq. (Q). As in Fig. [lj the left panel represents the EdS 
cosmology, the right panel the FlatLow cosmology. 



is chosen higher in the low-density universe, so that the pro- 
jected, physical distances of the measurement locations with re- 
spect to the cluster centre shrink. Closer to the cluster centre, 
the magnetic field is generally higher, and therefore these less 
magnetised clusters well reproduce the measurements. 

The overall structure of the magnetic field in the outer regions 
of the clusters differs from that near the cluster centre. Accre- 
tion and merger events arrange the magnetic fields in the outer 
cluster regions in coherent patterns on fairly large scales. Near 
the cluster centre, the gas flow pattern is almost randomised, 
and the magnetic field is consequently tangled and bent on fairly 
small scales. Although the magnetic pressure can be considered 
isotropic at or near the cluster centre, it may well be anisotropic 
further out. Even though the fraction of the magnetic relative 
to the thermal pressure, averaged within spherical shells, is ap- 
proximately constant across the cluster, the detailed effects of 
magnetic pressure on the gas flow depend on the direction of 
the magnetic field relative to other locally preferred directions, 
like the orientation of filaments surrounding a cluster, the local 
path of matter infall, the orbit of a merging sub-clump and such. 
Therefore, the same overall magnetic pressure can have differ- 
ent effects on the local dynamics of gas flows, depending on 
whether the magnetic field is ordered on scales comparable to 
the cluster scale, and on the orientation of the magnetic field. In 
addition, the magnetic field alters the dynamical time scale of a 
(simulated) cluster. Therefore, when comparing simulated clus- 
ters with and without magnetic fields at equal redshifts, we must 
be aware that we may be comparing them at slightly different 
evolutionary stages according to their intrinsic timescale. 

Figure || shows the temperature profile averaged across the 
set of ten simulated clusters. On the whole, the temperature in- 
creases towards the cluster centre by a factor of « 2. For dif- 
ferent cosmologies, the final clusters have a different amount 
of substructure. This leads to different shapes of the temper- 



ature profiles for the clusters in the two different cosmologies 
used. Although the detailed temperature profile of an individual 
cluster depends on its dynamical state and formation history, the 
tendency is clearly visible that the cluster temperature decreases 
with increasing initial magnetic field strength towards the cluster 
centre. Deviations from this general trend are possible in indi- 
vidual clusters because of the reasons mentioned before. The 
mean across the entire cluster sample, however, reveals a mono- 
tonic trend increasing towards the cluster centre. This behaviour 
is expected in view of the previous discussion about the impact 
of the magnetic field in different parts of the cluster. 

3.2. Temperature-mass relation 



Table 3. — This table shows the temperature difference be- 
tween the magnetised and non-magnetised clusters (columns 2 
and 4), and the rms (columns 3 and 5) of this distribution, for the 
three different initial magnetic field strengths. Columns 2 and 3 
are for the EdS, columns 4 and 5 for the FlatLow cosmologies, 
respectively. The temperatures are emission-weighted tempera- 
tures within the virial radius r2oo- All clusters shown in Figure [j] 
are used. 
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low 


4.3% 


2.7% 


2.0% 


2.1% 


medium 


3.1% 


3.7% 


2.4% 


2.1% 


high 


5.8% 


5.1% 


2.5% 


1.7% 



The effect of magnetic fields on the temperature-mass relation 
in galaxy clusters is shown in Fig. |[ Different initial magnetic 
field strengths are marked with different symbols. As shown be- 
fore, the effect of the magnetic pressure on the temperature is on 
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FIG. 3. — The large panels show the temperature-mass relation for the simulated cluster sample. Here, the temperature is the emission-weighted 
temperature within the radius ^200- All clusters are simulated with different initial magnetic field strengths marked by different symbols, as 
indicated in the plot. The right-hand overlays show the temperature change relative to non-magnetised clusters for the different initial magnetic 
field strengths in percent for each simulated cluster. The large scatter is obvious, but there is a general tendency towards lower temperatures in 
magnetised clusters. To follow this relation towards lower temperatures, smaller objects (marked as gray symbols) found in the simulations are 
also used, which lead to a total of 20 different clusters for the EdS cosmology (left panel) and 18 different clusters for the FlatLow cosmology 
(right panel). As low mass objects are represented by less particles, the resolution in these objects drops towards lower masses and the scatter in 
these objects is expected to be larger. 



average strongest near cluster centres. Therefore, the strength of 
the effect on the temperature-mass relation depends on the cho- 
sen radius within which the mean temperature is calculated. It 
changes also if we choose the emission-weighted instead of the 
mass-weighted temperature. For Fig. ||, we have chosen to cal- 
culate the emission-weighted temperature within the virial ra- 
dius r2oo in order to be comparable to the measurements. We 
also translated the masses of the simulated galaxy clusters to 
velocity dispersions. Figure [3] also shows the measured rela- 
tionship between cluster temperatures and velocity dispersion as 
dashed-dotted line. The fairly large scatter in individual clusters 
is clearly seen. In the mean, the temperature difference between 
the magnetised and non-magnetised clusters is 5%, and it can 
go up to 10 — 15% in individual clusters. 

Although the systematic effect of magnetic fields on the 
emission-weighted temperature profiles is small, the magnetic 
fields also increase the scatter in the mass-temperature relation. 
Values for the rms scatter are given in Table This effect is 
more pronounced in the EdS cosmology because of the smaller 
gas fraction in a high-density universe. 

4. Discussion 

We tested the influence of the non-thermal pressure support due 
to magnetic fields in galaxy clusters. To do so, we performed 
cosmological MHD simulations of ten clusters in two different 
cosmologies, each with four different initial magnetic field con- 
figurations, thus yielding a total of 90 cluster models. 

We chose initial magnetic field strengths such that the fi- 
nal intracluster fields were able to statistically reproduce the 
amplitude and spatial distribution of observed Faraday-rotation 



measurements. For such fields, we found a relatively small 
non-thermal pressure support of about 5%. The temperature 
change in the cluster cores is therefore on average w 5%, with 
the general tendency towards decreasing temperature with in- 
creasing magnetic field strength. This tendency is most pro- 
nounced in cosmologies with a low gas fraction, i.e. a high dark- 
matter density. The effect of the magnetic fields can be substan- 
tially stronger in individual clusters, introducing a cosmology- 
dependent scatter of up to ( 10 — 15) % into the mass-temperature 
relation. 

Our overall conclusion is that the statistical rotation-measure 
data, taken in the most straightforward interpretation, suggests 
that magnetic fields are generally dynamically unimportant in 
galaxy clusters. We are aware that simulations like those pre- 
sented here do not include the full range of physical processes 
required for a definitive calculation of intracluster fields, like 
small-scale turbulence and non-ideal magneto-hydrodynamics. 
A detailed study of such processes is currently out of reach for 
realistic numerical simulations, rendering our simulations illus- 
trative rather than definitive. We point out, however, that the 
inclusion of magnetic fields in cluster simulations does consti- 
tute an improvement over previous studies, even if it confirms 
earlier conjectures. 

5. ACKNOWLEDGEMENTS 

We wish to thank Harald Lesch for useful discussions, and the 
referee, Jean Eilek, for her very detailed and thoughtful report. 

REFERENCES 



Brio, M, Wo, C.C., 1988, J. Comput. Phys., 75, 400 

Clarke, T.E., Rronberg, P.P., Bohringer H., 1999, Proceedings of "Diffuse Ther- 
mal and Relativistic Plasma in Galaxy Clusters", eds. H. Bohringer, L. Fer- 
etti, P. Schuecker, MPE Report 271, 82, 1999 
Dolag, K., Bartelmann, M., Lesch, H., 1999, A&A, 348, 351 
Dolag, K., 2000, Proceedings of "Constructing the Universe with Clusters 

of Galaxie s^ TAP 7000 prk Flnrpnpp T~)iirrpf Rr Flnnipl Tiprhal aStTO- 

ph/0010149 . . 

Dolag, K., Ensslin, T., 2000, A&A in press, astro-ph/00 10296 
Kim, K.T., Kronberg, P.P., Tribble, PC, 199 1, ApJ. 379, 80 
Kronberg, P.P., 1994, Rep. Prog. Phys., 57, 325 
Lucy, L., 1977, AJ, 82, 1013 
Monaghan, J.J., 1992, ARA&A, 30, 543 
Steinmetz, M., 1996, MNRAS, 278, 1005 

Sugimoto, D., Chikada, Y., Makino, J., Ito, T., Ebisuzaki, T., Umemura, M., 

1990, Nat, 345, 33 
Vallee, J.P, MacLeod, J.M., Broten, N.W., 1986, A&A, 156, 386 
Vallee, J.P, MacLeod, J.M., Broten, N.W., 1987, ApL, 25, 181 
Wu, X., Fang, L., Xu, W., 1998, A&A, 338, 813 



